Velocity difference statistics in turbulence 
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We unify two approaches that have been taken to explain the non-Gaussian probability distribu- 
tion functions (PDFs) obtained in measurements of longitudinal velocity differences in turbulence, 
and we apply our approach to Couette- Taylor turbulence data. The first approach we consider 
was developed by Castaing and coworkers, who obtained the non-Gaussian velocity difference PDF 
from a superposition of Gaussian distributions for subsystems that have a particular energy dissi- 
pation rate at a fixed length scale [Castaing et al., Physica D 46, 177 (1990)]. Another approach 
was proposed by Beck and Cohen, who showed that the observed PDFs can be obtained from a 
superposition of Gaussian velocity difference PDFs in subsystems conditioned on the value of an 
intensive variable (inverse "effective temperature") in each subsystem [Beck and Cohen, Physica A 
322, 267 (2003)]. The intensive variable was defined for subsystems assuming local thermodynamic 
equilibrium, but no method was proposed for determining the size of a subsystem. We show that 
the Castaing and Beck-Cohen methods are related, and we present a way to determine subsystem 
size in the Beck-Cohen method. The application of our approach to Couette- Taylor turbulence 
(Reynolds number 540 000) yields a log-normal distribution of the intensive parameter, and the 
resultant velocity difference PDF agrees well the observed non-Gaussian velocity difference PDFs. 

PACS numbers: 47.27.Ak, 47.27.-i, 05.20.-y, 47.27.Jv 
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I. INTRODUCTION 

In Kolmogorov's 1941 theory (K41), the energy in fully 
developed three-dimensional turbulence cascades from 
large scales to small scales where it is dissipated pj. 
Turbulence in the cascade (the inertial range) is char- 
acterized by the probability distribution function (PDF) 
P(5v r ) for longitudinal velocity differences over a dis- 
tance r, 5v r (x) = e r ■ [v(x + r) — v(x)]), where e r is the 
direction of separation |2(. For r approaching the inte- 
gral scale where energy is injected, the PDF is Gaussian, 
while in the inertial range extending down to the dissi- 
pation scale ?7, intermittent large fluctuations lead to a 
non-Gaussian PDF with approximately exponential tails 

s 

Kolmogorov assumed a constant energy dissipation 
rate per unit volume, £ |. In 1944 Landau 3 
suggested that fluctuations of e averaged at scale r, 

E r (x,t) (= Jg +r e(x' ,t)dx'^j play a key role in turbulence. 
Such fluctuations were subsequently observed in many 
experiments [H 0, 0- I n 1962 Kolmogorov || and 
Obukhov 9] proposed a log-normal distribution of e r 
in the inertial range. The log-normal distribution was 
obtained in subsequent experiments and numerical sim- 
ulations e r [Hill mill. The non-Gaussian PDF of 
Sv r and the log-normal PDF of e r characterize turbulent 
flows. 

Different approaches have been taken by Castaing et 
al. 0] and by Beck and Cohen 01 to understand the 
non-Gaussian P(Sv r ). Castaing et al. assumed that sub- 



systems have different values of e r , but the subsystems 
have Gaussian PDFs of 5v r ; this assumption is supported 
by experiments. 

Beck and Cohen took a statistical mechanics approach, 
assuming that subsystems have a well-defined "effective 
temperature" , which for turbulent flow is identified with 
the variance of 5v r . The resultant P{5v r ) depends on 
the statistics of the distribution for the inverse effective 
temperatures in the subsystems. This dependence of the 
statistical distribution P(Sv r ) on the statistical distri- 
bution of subsystems led Beck and Cohen to call their 
approach super statistics |l5j . 

In this paper we note that the approaches of Castaing 
et al. and Beck and Cohen are both based on Bayes' 
theorem, 



P(x) = / P(x\y)P(y)dy, 
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(1) 



which is used to obtain the non-Gaussian P(Sv r ) from 
a conditional mixing of Gaussian PDFs in subsystems. 
However, the subsystems are chosen differently in the 
two approaches. 

We propose a method that does not require a determi- 
nation of e r from experimental data, nor does it require 
a fitting parameter to obtain the effective temperature 
PDF. We show that subsystems with Gaussian statistics 
can be chosen by examining moments of velocity differ- 
ence distributions in the subsystems. Our method, which 
involves no fitting parameters, leads to predictions for 
the non-Gaussian P(5v r ) that are in accord with data 
for turbulent Couette- Taylor flow [ljj . 

In Section [H] we present the Castaing et al. and 
Beck and Cohen methods, and in Section ITTT1 we describe 
the Couette- Taylor experiments and present results for 



2 



P(Sv r ). Section lIVl shows how subsystems can be sys- 
tematically chosen to obtain a prediction for P(Sv r ). The 
conclusions are presented in Section Ivl 



II. THEORY 
A. Method of Castaing et al. 

Castaing et al. 0] started with the observation from 
their experiments that that velocity difference distribu- 
tions for a given e r are Gaussian, and that e r is described 
by a log-normal distribution 13 13 13 HJ SH • The 
log-normal distribution for e r has also been obtained 
for e r in other experiments on fully developed turbu- 
lence [T^.l22 l l23.l2 ^ |. and in analyses of images of cloud 
patterns |25j . effective temperature fields in turbulence 
|2fij | , and magnetic fields in solar winds [2]} ■ 

To describe the evolution of P(Sv r ) from Gaussian at 
large scales to non-Gaussian at small scales I29ll30| . 
Castaing et al. proposed 0, IrH ] 



P{Sv r 



P(e r )P(5v r \e r )de r 



(2) 



The conditional PDF P(Sv r \s r ) in Eq. J5J is assumed to 
be a Gaussian distribution, P(5v r \e r ) = e - ^^ 2 '/( re r) 2li ; 
in accord with experimental observations [SHl l32l l33| . 
Kolmogorov @, Obukhov and Castaing 0] assumed 
a log-normal distribution of e r , 



P(e r ) 



X e (2n) 1 / 2 e r 



exp 



(me r - m s )' 
2Al 



(3) 



where m E and A e are respectively the mean and the stan- 
dard deviation of lne r . 

A difficulty in applying the approach of Castaing et 
al. is that energy dissipation rate at length scale r, e r , 
is not directly measured in experiments. By assuming 
homogeneous and isotropic conditions, e r (x) is defined 
as Ibv r (dv/dx) 2 dx. In practice, e r is determined 
from time series data, 



1 r N - 1 

1=1 



(Ax 



[v(xi+i) - v(xi)f 



(4) 



where Ax(= x 2 — xi) is the sampling separation the sum- 
mation i is over subsystems and xn — x\ = r 

SHU!. 

Even with this assumption, determination of e r (x) is 
difficult because of errors in evaluating the derivative 
from velocity data. Further error arises from the ap- 
plication of the Taylor frozen hypothesis at high frequen- 
cies IllllilillMlliliSlilllilillil. 



B. Superstatistics of Beck and Cohen 

Beck and Cohen's statistical approach considers a 
system far from thermodynamic equilibrium to consist 



of subsystems in local thermodynamic equilibrium |l5| . 
Each subsystem has a well-defined "effective tempera- 
ture" , but the subsystem effective temperatures need not 
be the same since the whole system is not in equilibrium. 
Beck and Cohen identify (5v r ) 2 with the kinetic energy 
of eddies of size r, E(Sv r ) — ^(5v r ) 2 , and the variance 
of dur is identified with an inverse effective temperature 
(3 [Ifj , given for a subsystem of size d by 



0d = 



((Sv r r) d -((6v r ) d )i-' 



(5) 



where (-)^ is an average over the size d. Then we have 

P(6v r ) = / P{(3 d )P{5v r \p d )d(3 d . (6) 
Jo 

where P(f3 d ) is the distribution of inverse effective tem- 
perature in subsystems of size d. 

A particular choice of P([3 d ), the \ 2 distribution, 
leads to the distribution associated with the nonexten- 
sive statistical mechanics of Tsallis, P(E) = (1 + (3(q — 
q is a parameter characterizing the 
nonextensivity [S(l+2) = S {iU_S£2) + {1 - q) S (1) ■ S {2) , 
where S is entropy function.] |l5ll46l|. A phenomenology 
similar to Beck and Cohen's was used in earlier oceano- 
graphic analysis that described the global non-Gaussian 
distribution of ocean surface velocity as a mixture of lo- 
cal Gaussians with x 2 -distributed variance 0, El] • The 
method of Beck and Cohen has been applied to fully de- 
veloped turbulence p^. IHol ] by introducing a fitting pa- 
rameter to determine the PDF of inverse effective tem- 
perature, rather than by directly measuring the PDF of 
inverse effective temperature. 

The Beck-Cohen method requires that the size d should 
be large compared to the distance r separating two 
points, and d should also be large enough so the subsys- 
tems contain enough data points to yield good statistics, 
but d must also be small enough so that subsystems are 
each described by a Gaussian distribution. Beck deter- 
mined the size of d using a fitting parameter involving 
the kurtosis of P(Sv r ) [5(1 ]. 



C. Unified view of PDFs 

The Castaing and Beck-Cohen methods are similar ex- 
cept in the way they divide a system into subsystems. 
Castaing et al. sample velocity differences conditioned 
by the averaged energy dissipation rate e r , while Beck 
and Cohen use velocity differences conditioned by the in- 
verse effective temperature j3 d - Castaing et al. need one 
fixed length scale, the separation distance r between two 
points; 5v r and e r are defined at this scale and are re- 
lated through Bayes' theorem. The Beck-Cohen method 
involves two length scales, the distance r separating two 
points and the size d of the subsystems in the statistical 
analysis. 
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The Castaing and Beck-Cohen methods can be con- 
nected if the two conditioning variables (e r and (3 d ) are 
correlated. Using Eq. (J2J and Bayes' theorem, we con- 
vert Castaing's method into Beck-Cohen's method, 



P(Sv r ) 



P(5v r \e r )P(e r )de r 



oo poo 



P(Sv r \l3 d )P(/3 d \s r )df3 d P(e r )de r 



JO 

P(Sv r \f3 d ) 



P((3 d \e r )P(e r )de r 



P(6v r \p d )P(/3 d )dp d . 



d/3 d (7) 
(8) 



Now, let's assume a log- normal distribution of (3 d at the 
fixed e r , 



P{(3 d \e r ) oc -j-exp 
Pd 



(ln/3(j — ah\e r Y 



2X t 



(9) 



where A t is the standard deviation of In f3 d conditioned 
to e r , and a is a parameter. Using Eqs. ©, 0, and 10, 
we have 

(In (3 d - a\ne r f 



P{5v r ) oc / P(5v r \f3 d ) I exp 
Jo Jo 

( (lne r — m £ ) 2 

xexp — 



2V 

d(\n e r )d(\n(3 d ) 



oc 



P(5v r \p d )exp ( 2 m)2 ) dQnfa) 



(10) 



Thus with the assumption of a log- normal distribution of 
(i d conditioned on e r , we have that Castaing's method is 



equivalent to Beck-Cohen's method. In Section llV Dl the 
log-normal PDF of P((3 d \e r ) is verified in experiments. 



III. EXPERIMENT 

We describe here an experiment on turbulent Couette- 
Taylor flow by Lewis and Swinney |16l l51| , and in the 
next section we will analyze data from this experiment 
to deduce P{(3) and a prediction for P(5v r ). The fluid 
was contained in the annular region between two concen- 
tric cylinders with an inner radius of b — 22.085 cm and 
an outer radius of a = 15.999 cm; thus the ratio of in- 
ner to outer radius was 0.724. The height of the annulus 
was 69.5 cm, which yields a value of 11.4 for the ratio of 
height to the gap. The inner cylinder angular rotation 
rate fi was 8 x 2tt rad/s; the outer cylinder was at rest. 
The ends of the annulus rotated at the same rate as the 
inner cylinder. The fluid was water with a viscosity v 
of 0.00968 cm 2 /s at the working effective temperature. 
Defining the Reynolds number as Re = Qa(b — a)/v 
yields for the Reynolds number 540 000 fl6|. 

A hot film probe was used to measure the time depen- 
dence of the azimuthal component of the velocity in the 




FIG. 1: An example of the Couette- Taylor velocity difference 
data, obtained by subtracting velocities at two points with a 
separation r — 46r; = 0.134 cm, where r\ is the Kolmogorov 
length scale. The inset shows the velocity differences on a 
finer length scale. 



center of the gap at a distance 4.35 cm above mid-height 
of the annulus. The Taylor frozen turbulence hypothe- 
sis was used to convert the velocity time series data to 
velocity field data. The turbulent intensity (the ratio of 
the root mean squared velocity to the mean velocity) was 
less than 6%. 

The uncertainties shown on our graphs correspond to 
the standard deviation of 20 independent experiments. 
The velocity measurements were made with a sampling 
rate 2500 times the inner cylinder rotation frequency; this 
corresponds to a spatial separation of 0.017 cm between 
successive velocity values. The longitudinal velocity dif- 
ferences Sv r that we analyze are for points separated by 
a small distance, r = 0.134 cm, where the probability 
distribution function has approximately exponential tails 
|l6j | . An example of the measurements of Sv r (t) is shown 
in Fig. ^ The separation r = 0.134 cm corresponds to 
46?y, where rj is Kolmogorov scale [3. (TheKolmo gorov 
dissipation scale was obtained by calculating the dissipa- 
tion from energy spectra: rj = (i^/e) 1 / 4 , where the dis- 
sipation rate is given by e — \hv J k 2 E(k)dk [Rj-) The 
window size d we use for determining the local inverse ef- 
fective temperature (3 is typically 0.9 cm, nearly an order 
of magnitude larger than the value of r. 



IV. RESULTS 

A. Probability density function of inverse effective 
temperature 

Several distributions for inverse effective temperature 
(3 d have been discussed by Beck and Cohen . Here we 
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FIG. 2: Comparison of \ 2 an( A log-normal distributions to the experimental distribution for inverse effective temperature in 
subsystems of size (a) rf=0.9 cm and (c) d=3 cm. The dash-dotted lines represent the \ 2 distribution, and the solid line 
represents the log-normal distribution; both have the same mean and variance as the 20 independent experiments (error bars 
correspond to one standard deviation). The panels on the right, (b) and (d), show the difference between the experimental 
PDF for (3d and the \ 2 (plus signs) and log-normal (bullets) distributions for (a) d=0.9 cm and (c) d=3 cm. The shaded area 
represents the experimental uncertainty (standard deviation of 20 experiments). 



consider the log normal and y 2 distributions, which are 
most applicable to turbulent flow. Due to multiplicative 
processes in turbulence, the log-normal distribution is 
often observed for positive-definite quantities (such as e r ) 
El 13 13 13 IU H3 A log-normally distributed (5 d is 
given by 

p{Pd) = ^ww d exp { — ^ — J 



where s = y ln(l + cr^ d /Pd ) an( l m = 

\og((3d 2 / \j Pd 2 + & 2 p d ) are parameters, and fid and 
<j/3 d are respectively the mean and standard deviation of 
Pd- 

The \ 2 distribution of /3d is given by 
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FIG. 3: The parameters %- (circles) and m (triangles), ob- 
tained from fits of the inverse effective temperature (3d (de- 
duced from Couette- Taylor turbulence data) to a log-normal 
distribution, as a function of subsystem size d (see Eq. I1H . 
(s and m are the variance and mean of logarithmic inverse 
effective temperature.) The parameters %- and m are ap- 
proximately equal (see text) and are described by a power 
law, m oc d~ 3//4 (solid line). 




d (cm) 



FIG. 4: The dependence of the third and fourth moments of 
P(Sv r ) on the size d of the subsystems. For sufficiently small 
d, P(5v r ) should be Gaussian, which means the values of the 
third and fourth moments should have the values zero and 
three, respectively. We find that at d « 0.9cm, the conditional 
distribution Eq. (111! is close to Gaussian (see text). 



where c — [3d / a p d an d b — er 2 ^ / (3d and V is the gamma 
function. The \ 2 distribution has been observed in re- 
cent measurement of wind turbulence |5^ | . The statisti- 
cal properties of different distributions are discussed in 

The experimental PDF for (3d is compared in Fig. [3 
with a log-normal distribution and with a x 2 distribution 
for two subsystem sizes d, 0.9 cm and 3 cm. The mean 
(3d and variance cr| d of the inverse effective temperature 
determine the parameters s, m, b and c. For small d, the 
log-normal and % 2 differ significantly, but for large d they 
become closer together [Fig. |2Jc) and (d)]. The decrease 
in variance of (3d with increasing d is similar to decrease 
observed in the variance of e r with increasing r |53| . 

The difference between the PDF of (3d from experiment 
and the \ 2 an d log-normal distributions is shown in Fig. 
|2Jb) and (d). For d=0.9 cm, the log-normal distribution 
fits the data within the experimental uncertainty except 
small (3d regions, while the x 2 distribution deviates from 
the observations by an amount that is large compared to 
the uncertainty. For d=3 cm, the log-normal distribution 
fits the distribution of (3d whereas the x 2 distribution 
does not. 

The log normal distribution involves two parame- 
ters, s and m, which depend on subsystem size, as shown 

in Fig. |21 This figure suggests a relationship between s 

2 

and m, m — 4^-, which is supported by a calculation in 
Castaing et al. (see Section 4.3.1 in jT4|). 



B. Conditional probability and the proper 
subsystem size 



In the statistical approach of Beck and Cohen, the sub- 
system size d should be sufficiently small so that P((3d) is 
Gaussian, corresponding to local thermodynamic equilib- 
rium in the subsystems. However, in practice the d — > 
limit is inaccessible because as d becomes very small, the 
number of data points becomes too small to allow ac- 
curate determination of the variance of (3d- So what is 
optimal choice of dl We address this question by exam- 
ining the third moment (skewness) and fourth moment 
(kurtosis) of Sv r , which should be equal respectively to 
zero and three for a Gaussian distribution. In principle 
we could also examine fifth and higher moments, but be- 
cause of the sensitivity of the higher moments to noise, we 
limit our considerations to the third and fourth moments. 
Plotting the third and fourth moments as a function of 
d, as shown in Fig. ^ we find that the optimal value of 
d for our data is 1.0-1.2 cm, which is the only range in 
which the kurtosis is approximately given by the value 
for a Gaussian. The skewness is small and negative for 
d > 0.5 cm, but becomes strongly positive for d < 0.5 cm, 
reflecting a cascade of energy to smaller length scales. We 
conclude that d=0.9 cm is the optimal subsystem size for 
our data. 
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P(8v. 




Sv r /(8Vr 




Sv r / (6v r )i 



FIG. 5: Comparison experimental results (dots) for P(8v r ) 
with the prediction of the Beck-Cohen method for a subsys- 
tem with the optimal size of 0.9 cm (bold line) on semi-log 
scale in (a) and relative error between theoretical and exper- 
imental values in (b). For comparison, we also show in (a) 
and (b) the predictions for subsystems of size 0.3 cm (thin 
dashed line) and 3 cm (thin dash-dot line) and a Gaussian 
distribution (dashed line). 



lne r 



4 5 

In /3 d 



FIG. 6: The relation between (id and e r . The solid vertical 
lines represent standard deviations at a fixed fid and the dots 
represent the mean values. The dashed line is fid oc (£,-)~ 2//3 . 



P(5v r ) 



1 

2ns 
x exp 



d(3 d f3 d 1/2 exp 



(log/3rf - m) 
2s 2 



(13) 



where s and m are determined from experiment for the 
optimal subsystem size d. There is no explicit form for 
the improper integral in Eq. 1(13(1 so we evaluate the inte- 
gral numerically, using the limits ([min/3,2, max/J^]) mea- 
sured in experiments instead of the theoretical integral 
domain, [0, oo). 

The results for P(5v r ) obtained by numerical integra- 
tion of Ijl 8(1 are shown in Fig. The data are described 
much better by the predicted probability distribution 
than by a Gaussian. The observed approximate power 
law tails are similar to the predicted distribution func- 
tion. 



D. Castaing and Beck-Cohen methods 



C. Probability distribution of 8v r 



We found a log-normal distribution of fid fits the tur- 
bulence data over a wide range in d (Sect ion llV A (I . With 
the log normal distribution of (3d for the optimal value 
of d (0.9 cm, Fig. and the conditional Gaussian dis- 
tribution of 5v r for that (3d, we obtain the probability 
distribution of Sv r by the method of Beck and Cohen, 



If the two conditioning quantities in the Castaing and 
Beck-Cohen methods (e r and (3d, respectively) are cor- 
related as a power-law, through Bayes' theorem the two 
methods can be seen to be the same (see Eq. (|5J|). With 
the surrogate definition of e r as in Eq. and a proper 
subsystem size (Section IIVB(I . we find that (3d and e r 
exhibit a power-law relation, as Fig. illustrates. In this 
sense, the Castaing and Beck-Cohen methods describe 
the same PDF of 8v r through the different conditional 
values which are correlated. Our experimental observa- 
tion of a relation (3d oc (e r )~ 2 / 3 in Fig. follows also 
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ln/V(ln/3 d ) r . 



where square brackets [•] denote the dimension of a phys- 
ical quantity, T is the dimension of time and L is the 
dimension of length. 

The probability of {3d conditioned to e r , P(j3d\£ r ), is 
log-normally distributed, as Fig. [7\ illustrates. Our as- 
sumption in Eq. J^J holds with the surrogate e r and (3d, 
where d is properly chosen (Section llV fj|) . Thus the inte- 
gral of two log- normal distributions, J P((3d\£ r )P{£r)de r , 
is another log- normal distribution, P((3d). That is, if 
P{(3d\£r) is a log-normal distribution with the mean 
of lne r , a log-normal distribution of e r in Castaing's 
method is equivalent with a log-normal distribution of 
fld in Beck-Cohen's method. 




-5 5 



ln/3 d /(ln/3d)r 



FIG. 7: The Gaussian distribution of In (3d conditioned by 
e r , plotted on (a) log and (b) linear scales. The solid lines 
represent a Gaussian distribution of In (3d, that is, the log- 
normal distribution of /3d. The dots represent the mean values 
of In (3d from experiments. 



V. CONCLUSIONS 



Both Castaing and Beck-Cohen methods have been 
very successful in describing the non-Gaussian distribu- 
tion of velocity differences in turbulence 0, ISfj. Al- 
though the relation of Beck-Cohen's method and Tsallis 
statistics |6lJ to turbulence has been questione d loOl lodl 
Ulil, the fit to data is quite good [H |H |U iof We 
have presented a method for determining subsystem size 
in the Beck-Cohen method, thus eliminating the need for 
a fitting parameter. 

We have also shown that Castaing's method can be 
converted to Beck-Cohen method - the log-normal dis- 
tribution of e r in Castaing's method gives rise to a log- 
normal distribution of (3d in Beck-Cohen's method. In 
that sense, the two methods describe the non-Gaussian 
distribution of Sv r in the same way, P(5v r ) = J Gaussian 
distribution x log-normal distribution. 



from a dimensional analysis, 



=>0 d oc r- 2 / 3 e- 2/3 , (14) 
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